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During the last few years the Tevatron has improved the bounds on 
rare B-meson decays into two leptons. Sensitivity to this decay is also one 
of the benchmark goals for LHCb performance. We compute the complete 
1-loop MSSM contribution to B® d — > and study the predictions for 

arbitrary flavour mixing parameters. We discuss the possibility of both 
enhancing and suppressing the branching ratios relative to their SM expec- 
tations. We find that there are "cancellation regions" in parameter space 
where the branching ratio is suppressed well below the SM expectation, 
making it effectively invisible to the LHC. 

PACS numbers: 12.60.Jv; 12.15.Ff; 13.20.He 

1. Introduction 

One of the most promising signals for new physics at the LHC is the rare 
decay B® — * /i + fi~ . The decay is strongly suppressed but especially 'clean' 
because its final state is easily tagged and its only hadronic uncertainties 
come from the hadronic decay constant fg g [T]. The LHC will be the first 
experiment to be able to probe this decay channel down to its Standard 
Model (SM) branching ratio. In particular the LHCb will be able to directly 
probe the SM predictions at 3a (5a) significance with 2 fb _1 (6 fb _1 ) of 
data, or after about 1 year (3 years) of design luminosity [2]. It is not clear 
whether LHC can reach the SM expectation for B d — > n + n~ ■ 

The current experimental status and the SM predictions for the branch- 
ing ratios B® d — ► are: 



Channel 


Expt. 


Bound (90% CL) 


SM Prediction 


B" s -> fj+fT 
B° d - 


CDF II [3] 
CDF II [3J 


< 4.7 x 10~ 8 

< 1.5 x 1(T 8 


(4.8 ± 1.3) x l(T y 
(1.4 ±0.4) x 10" 10 



At the dawn of the LHC era, it is important to understand the possible 
contributions of new physics to the B® d — > decays. They could be 



* Presented at the FlaviaNet Topical Workshop, 23-27 July 2009, Kazimierz, Poland. 



(1) 



2 



flavianet'jr printed on November 17, 2009 



particularly large in the Minimal Super symmetric Standard Model (MSSM). 
Under the assumption of large tan (5 and Minimal Flavour Violation (MFV, 
flavour violation given by CKM matrix only), the branching ratio for B® — > 
is dominated by the Higgs penguin mode and can be significantly 
enhanced over the SM expectation, as can be seen from the approximate 
formulae (for detailed study focusing on the resummation of tan/3 see [1]): 

With the upcoming new experimental probes of B(B® d — > it is 

important to perform a full calculation of this decay rate without a priori 
assumptions on the pattern of flavour mixing. In particular, in the region 
of low tan (3, the interference of box, Z- and Higgs-penguin diagrams could 
conceivably lead to both enhancement (even visible at Tevatron) or a can- 
cellation that would suppress the branching ratio below the SM prediction 
(making it invisible at the LHC). The status of this decay could become an 
important factor for planned LHCb upgrades and determining whether they 
should focus on increasing sensitivity to B® — > fJ, + fJ,~ or instead reach for 
the smaller branching ratio of B® — > [i + fi~. The relevant calculations have 
been presented in details in [5], we stress here the most important results. 



2. Effective Operators and Branching Ratios 

The effective Hamiltonian for the quarks-to-leptons transition q I q J — > 
£ +K £~ L , with q 1 = d,q 2 = s,q 3 = b and I 1 = e, £ 2 = fi, £ 3 = r, reads: 

T~L = , , ( CvxyOvxy + CsxyOsxy + CtxOtx), (2) 

where flavour and colour indices have been suppressed for brevity. The 
(V)ector and (S)calar operators are respectively given by 

Osx K y L = (jpPx^WPY^). (3) 
The explicit form of the Wilson coefficients in the MSSM is given in [5]. 
The (T)ensor operator contributions and photon penguin contribution to 
Bg — > l + £'~ vanishes in matrix element calculations. We do not consider 
the very large tan [3 > 30 scenario, thus no resummation of higher orders in 
tan/3 is necessary. The matrix element for B® d — > £ + £'~ decay is: 

M = F s ££ + F P l lb £ + F vV n 1(i £ + F Al H W l, (4) 

where the £s correspond to external lepton spinors. The (S)calar, (Pseudo- 
scalar, (V)ector and (A)xial-vector form factors in Eq. (JH) are given by 
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Fs = 7 (Csll + Cslr - Csrr - Csrl), (5) 

4 m b + m s{d) 

I M| d f Bs{d) 

Fp = (-Csll + Cslr- Csrr + Csrl), (6) 

4 m b + m s (d) 

Fy = --fB s(d) (CvLL + CvLR- Cvrr- Cvrl), (7) 

1 

Fa = — ^/B s ( d ) (-CVll + Cvlh - Cyfifl + Cvrl)- (8) 

The general expression for B(B® d — » is rather complicated [5]. 

For the most important — > fi + n~ decays it simplifies greatly and reads 
approximately as (q = s, d) 

- ^^(l^l 2 + |Fp + 2m,F A | 2 ), (9) 
where rs 9 is the lifetime of B q meson and we have taken the limit * 0. 

3. Numerical Analysis of S SjC j — > l^ + fJ-~ 

We may distinguish two possible scenarios for the relative size of the 
MSSM contributions to the right-hand side of Eq. (|9|): 

1. Higgs penguin domination or large tan /3 > 10. In this regime one can ex- 
pect large Fs, Fp and an enhancement of the branching ratios as in Eq. (pQ) 
(barring some possible GIM-type cancellations leading to i 7 fp SY ~ |6j). 

2. Comparable Box, Z- and Higgs-penguin contributions or low tan (5 < 10. 
In this case the super symmetric Higgs-mediated form factors Fs p may be- 
come comparable to or even smaller than Fa- Either an enhancement or a 
suppression of the branching ratios is possible depending on the particular 
choice of MSSM parameters. 

An enhancement of the branching ratios occurs generically in most of the 
MSSM parameter space. It is a bit trickier to suppress the branching ratios 
below their SM predictions. This is the case we would like to investigate 
further. We would like to find the minima of B(B® d — > fi + fi~), i.e. the 
minima of Eq. ©. We distinguish between two cases: 

a) F P + 2m e F A w and F P » F s , (10) 
a) \Fs\ ~ \Fp\ w \Fa\ ~ 0. (11) 

In the case a), the pseudoscalar and axial contributions cancel while 
the scalar contribution is negligible. The case b), happens when Higgs 
contributions are negligible compared to the axial contribution (i.e. low 
tan (3 and large Ma) and Fa becomes small due to cancellations among the 
Cvxy coefficients in Eq. ([8]). Our numerical analysis shows that both can- 
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Table 1 


Parameter 


Symbol 


Min 


Max 


Step 


Ratio of Higgs vevs 


tan (3 


2 


30 


varied 


CKM phase 


1 





7T 


7r/25 


CP-odd Higgs mass 


M A 


100 


500 


200 


SUSY Higgs mixing 




-450 


450 


300 


SU (2) gaugino mass 


M 2 


100 


500 


200 


Gluino mass 


M 3 


3M 2 


3M 2 





SUSY scale 


Msvsy 


500 


1000 


500 


Slepton Masses 




MguSY/3 


A^susy/3 





Left top squark mass 


Mr, 
Ql 


200 


500 


300 


Right bottom squark mass 


Mr 
or 


200 


500 


300 


Right top squark mass 


Mr 

t T1 


150 


300 


150 


Mass insertion 


$13 $23 
°dLL' dLL 


-1 


1 


1/10 


Mass insertion 


$13 $23 


-0.1 


0.1 


1/100 



cellations are possible but require a certain amount of fine tuning once 
constraints from other FCNC measurements are imposed. 

To quantitatively study the effects mentioned above we perform a scan 
over the MSSM parameters, not restricted to MFV scenario. Flavour vio- 
lation is parametrised by the "mass insertions" , defined as in [8] 



°QXY 



L Q)XY 



(12) 



where I, J are squark flavours, X, Y denote field chirality, and Q indicates 
the up or down squark sector. Note that we use the <5's for presentation only, 
not as expansion parameters - we numerically diagonalize all mass matrices. 

The ranges of variation over MSSM parameters are shown in Table 1, 
where "SUSY scale" denotes the common mass parameter for the first two 
squark generations, tan /3 takes on values (2, 4, 6, 8, 10, 13, 16, 19, 22, 25, 30), 
^dLL' ^dLR'A 4 '^ are taken to be real and the trilinear soft couplings are 
set to A t = A b = Mfi and A f = M 



L ' 



To realistically estimate the allowed range for B{B S 



one 



must account for the experimental constraints from other rare decays. For 
that, we use the library of numerical codes developed in the framework 
of the general MSSM in [1E1 [TOj and take into account the set of 
observables listed in Table ^J. In our scan, for the quantities in Table 2 we 
require (depending if the experimental result or only the upper bound is 
known) 

\Q ex P -Q th \< 3AQ ex P + q\Q th \, 
or (13) 
(l+q)\Q th \ <Q exp . 



1 For Higgs mass we use LEP data [11] > 92.8—114 GeV depending on sin 2 (a — (5). 
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Table 2 



Quantity 



Current Measurement Experimental Error 




> 46 GeV 

> 94 GeV 

> 89 GeV 

> 95.7 GeV 

> 92.8 GeV 



\zk\ 
I AM* | 
| AMo | 
AM Bd 
AM Bs 

Br(B —> X s j) 

Bt(K L -> 7T°^) 

Br(if+ -> tt+j/P) 
Electron EDM 
Neutron EDM 



2.232 • 10~ 3 0.007-10- 3 

3.483 • 10~ 15 0.006 • 10~ 15 

< 0.46- 10~ 13 

3.337 • 10~ 13 GeV 0.033 • 10~ 13 GeV 

116.96 ■ 10~ 13 GeV 0.79 ■ 10" 13 GeV 

3.34 -10" 4 0.38 -10~ 4 

< 1.5 • 10~ 10 

1.5 ■ lO" 10 1.3 ■ 10~ 10 



< 0.07- 10" 26 

< 0.63- 10~ 25 



3AQ exp and \qQ th \ in Eq. (|13)) represent the 3a experimental error and the 
theoretical error respectively. The latter differs from quantity to quantity 
and is usually smaller than the value q = 50% which we assume generically 
in all calculations. The increased "theoretical error" is used to account for 
the limited density of a numerical scan, simultaneously avoiding strong fine 
tuning between MSSM parameters (for detailed discussion see |10|). 

Fig. [1] shows the predictions for B(B® — > /x + /i~) over a general scan of 
20 million points in parameter space, including the bounds of Table 2. We 
vary S^ L (upper panel) and 5^ LR (lower panel) one at a time while setting 
the other to zero (results are also weakly sensitive to other S's). When 5^ L 
is varied in the range [—1,1], we find B(B® — > ) m j n ~ 10~ 9 . This 

minimum is almost independent of tan (5 but depends on the magnitude of 
l^dizJ- ^he u PP er bound set by CDF, depicted as a solid red line, can be 
attained even with very low values of tan j3. 

More interesting is the case when S^lr ^ s var i e d in the range [—0.1, 0.1]. 
We find a narrow cancellation region around sa —0.01 and tan (5 < 10 
where B(B® — > fJ. + fJ-~) m in ~ 10~ 12 (lower right panel). This is three orders 
of magnitude lower than the SM prediction, making it effectively unobserv- 
able at the LHC. In order to better understand cancellation region we study 
a representative point with a very low branching ratio (all masses in GeV): 

tan/3 = 4, Ma = 300, ^ = -450, M 2 = 100, M 3 = 300, 
SUSY scale = 400, M~ tR = 150, A tJb = M ih = M~ b ^ r = 600. (14) 

The 'Box', 'Higgs' and 'Z' lines in Fig.[2]indicate the value of B{B° S -> n + fi~) 
given by only the listed contribution with all others set to zero. In the 
cancellation region the box contribution is negligible while the Higgs- and Z- 
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Fig. 1. Upper panel: Predictions for B(B S — > versus tan/3 (left) and (^ll 

(right) from the scan of MSSM parameters. The upper solid line shows the current 
upper bound from the Tevatron and the lower dashed line the SM expectation. 
Lower panel: Similar to the upper panel but with S^lr var i e d- 

penguin magnitudes are comparable. Thus the latter should cancel, which 
we illustrate individually plotting the absolute values of the form factors 
Fs p and 2m fjL FA of Eqs. (J5HE|) in the right panel of Fig. [2j At the minimum 
point of the total branching ratio \Fp\ is approximately equal to |2?rtui^| 
and \Fg\ is negligibly small. This can be explained from the form of Eqs. ([5]) 
and (O. If one assumes <53? » = ($dLR) > then Cslr and Csrl, the two 



Wilson coefficients most sensitive to the variation of S^p, have similar sizes 
and opposite sign, interfering destructively in the amplitude. 

We performed similar scan also for the meson decay, B& — » 
varying 5^ LL or 5^ LR instead of 5 2 d \ L or S^ R along with the other SUSY 
parameters. For both cases there exist points where B(B ( i — > /i + /i~) is 
reduced by an order of magnitude relative to the SM. These points are more 
sensitive to low tan/3 in the 'LL' case and fall into the case of Eq. (jlip . On 
the opposite, the ratio R = B{Bd — ► /i + fj,~)/B(B s — > fi + n~), which in the 
SM is fixed to R « |W^| 2 < 0.03, in the MSSM can be enhanced by a 
factor of 10 even for small values of S\\ L or S^ R . This suggests that collider 
searches for B{Bd — > are as important as those for £>(i? s — > 

Our results lead also to bounds on 5 parameters. As can be seen from 
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-0.015-0.01 -0.005 0.005 0.01 0.015 -0.015 -0.01 -0.005 0.005 0.01 0.015 



Fig. 2. Contributions to B(B° d — » Ior the parameters in Eq. (fl4|) versus 

"dLR- Left: Contributions from the various diagram classes. Right: Magnitude of 
the form factors appearing in Eqs. ([SHH]). 

Figs [[] and scan results for B(Bd — > fi + fi~), SdLL^dLL are weakly con- 
strained, they can take on values up to ~ 0.9 and still pass the constraints 
in Table 2, though points beyond 0.3 are less dense. Bounds in the 'LR' 
sector are tighter, 5f LR ,8f LR < 0.08. 

4. Conclusions 

We have discussed results of a complete, 1-loop calculation of the branch- 
ing ratios for the rare decay modes B® , — > without resorting to the 
limits of large tan (3 and MFV scenario and performed a numerical explo- 
ration of the MSSM parameter space. We find that there exist cancellation 
regions where the contribution of diagrams with supersymmetric particles 
interferes destructively with purely SM diagrams, thus allowing the branch- 
ing ratio to be significantly smaller than the SM prediction. We identify 
possible mechanisms of such cancellations and explain why they can occur 
for certain regions of parameter space. Such effects may effectively hide 
the dimuon B® decay mode from the LHCb even though it is supposed to 
be one of the experiment's benchmark modes. We have also shown that, 
barring the cancellations mentioned above, supersymmetric contributions 
in the general MSSM typically tend to enhance the branching ratio for 
B® d ^> even for moderate values of tan/3 < 10 so that an experimen- 
tal measurement close to the SM prediction would put strong bounds on the 
size of allowed flavour violation in the squark sector. Finally, we show that 
the B^ — > decay can also be either suppressed or enhanced compared 
to its SM expectation, leading in some cases to a situation where the rate 
of the B® decay is larger then that of the B®. 
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